DOMAIN OF ATTRACTION OF GAUSSIAN 
PROBABILITY OPERATORS IN QUANTUM LIMIT 
; THEORY 

O 

^ ■ KATARZYNA LUBNAUER AND ANDRZEJ LUCZAK 

> 

O ■ Abstract. We characterise the class of probabiUty operators be- 

I longing to the domain of attraction of Gaussian limits in the setup 

• which is a slight generalisation of Urbanik's scheme of noncommu- 

I tative probability limit theorems. 

< 

^ I In our investigation of the domain of attraction in quantum limit 

theory we adopt the approach introduced in the fundamental paper [6J 
which can be briefly described as follows. Let 7i be a separable Hilbert 
space. By a probability operator we mean a positive operator on Ti. 
^ I of unit trace. It is well known that such operators are in a one-to-one 

correspondence with normal states p on B(7i), and this correspondence 
t::J- . is given by the formula 

p{A)=tTAT, AeM{n). 

^ . The set of all probability operators on H will be denoted by ^. By 

! £^ we shall denote the set of all trace-class operators on H, and by 

£^ — the set of all Hilbert-Schmidt operators. 



1. Preliminaries and notation 



Let z V{z) be an irreducible projective unitary representation of 
the group M^'^ on Ti, satisfying the Weyl-Segal commutation relations 

(1) V{z)V{z') = e-^^^''''^V{z + z'), 

where z, z' G z = {xx, yi, x^, y^), 2;' = {x\,y'^, ■■■,x'^, y'^), and 

d 

k=l 

Fix z G M^*^. It is easily seen that {V{tz) : t G M} is a one parameter 
unitary group, thus by the Stone theorem there is a selfadjoint operator 
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R{z) on H such that 

V{tz) = e''^'^'\ teR, 

consequently, 

V{z) = e'^^'l 

The operators R{z) are called in |3J canonical observables. Let 



/oo 
\E,{dX) 



be the spectral representation of R{z). For a probability operator T 
we define its mean value mj{z), second moment m^{z) and variance 
cr^{z) by the formulae 



/oo 
A ti TE,{dX) 
-oo 

/oo 
iiTE,{d\) 
-oo 



-oo 

/oo 
(A - mf (2))2 trTE^(dA) = ml{z) - ml{zf 
-oo 

(cf. [3], Chapter V, $ 4]). Note that the notions defined above corre- 
spond to the mean value (expectation), second moment and variance, 
respectively, of the Borel probability measure /iz determined by the 
formula 

(2) /i,(A) =trm(A), AgS(M). 

A probability operator T is said to have finite variance if for each 
z G M^'^, ct'^{z) < oo (equivalently, rn^i^z) < oo). 

For a probability operator T we define its characteristic function 
f: ^ C as 



(3) T{z) =tiTV{z), zE 



b2d 



T has the following property called A-positive definiteness: for arbi- 
trary complex numbers Ci, . . . , c„ and vectors zi, . . . , Zn E M^'^ 

n 

i,fc=i 

'Quantum Bochner's theorem' states that for a complex- valued func- 
tion /: M?'^ — > C we have f = T for a certain probability operator 
T if and only if / is A-positive definite, continuous at the origin and 
/(O) = 1 (cf. m Chapter V, $ 4]). 
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It is immediately seen that for an arbitrary probability operator T 
and an arbitrary Zq G M?'^ the function 

is the characteristic function of some probability operator. 

Formula ([3]) for T G defines a map which extends uniquely to 
a linear isometry from onto the space of all complex- valued square 
integrable with respect to Lebesgue measure functions / with the norm 




(cf. Chapter V, $ 3, Theorem 3.2]). 

Let 2t be the set of all Hilbert-Schmidt operators T for which T 
vanishes at infinity. We define the convolution * in 21 by setting 

Moreover, we put ||T|| = Then 

\\T^^T2\\ ^ IIT1IIIIT2II, 

and consequently, the convolution algebra 21 is a Banach algebra with- 
out unit. The following inclusions hold true 

<p c i:^ c 21 c £^ 

(cf. 0). 



2. Statement of the problem 

The general scheme of quantum limit theorems introduced in [6] is 
as follows. For a triangular array {Tkn '■ k = n = 1,2...} 

of probability operators, a norming array {a^n '■ k = 1, . . . , A;„; = 
1, 2, . . . } of positive numbers, and a sequence {zn} of elements from R^'^ 
we form probability operators Sn defined by the characteristic functions 

(4) Sn{z) = e'^'"^'^l[nn{aknz), z E M^". 

k=l 

The norming constants a^n should satisfy the assumption of admissi- 
bility which means that the maps 

R^"" 3 z^Y[fk{aknz) 

k=l 
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are the characteristic functions of some probabihty operators for each 
n and any probabihty operators Ti, . . . , r„. Now if 



hm S„{z] 



S{z) 



z e 



t>2d 



for some function S, then from quantum Bochner's theorem it foUows 
that 5* is the characteristic function of some uniquely determined prob- 
abihty operator S. In this case S is caUed the limit operator. In the 
paper [6] the class of limit operators was described under the assump- 
tion of uniform infinitesimality of the operators from 21 given by the 
functions {T^nicLkn-) : k = 1, kn', n = 1,2, .. .}, analogously to the 
case of the classical infinitely divisible limit laws, while in the paper [S] 
for the case d = 1 norming by arbitrary 2x2 matrices was considered. 
We shall be concerned with a quantum counterpart of the classical sta- 
ble limit laws, i.e. we assume that fc„ = n and Ti„ = ■ ■ ■ = T„„ = T 
for some probability operator T. As for norming we adopt the above- 
mentioned more general approach and as the norming matrices we take 
matrices An of the form 



(5) 



An) 










An) 











An) 
''d 










As in the scalar case we put the assumption of admissibility of the 
matrices A^ which means that for each n and any probability operators 
Ti, . . . ,Tn the function 



t,2d 



3 2; f-^ JJ fk{Anz) = Y[ fk{a'C^xi,a'C^yi, a^^^Xd, a^d^yd) 



k=l 



k=l 



is the characteristic function of some probability operator. In this case 
the limit operator S will be said to belong to the domain of attraction 
of the probability operator T. 

To justify this approach let us look at the fundamental notion of the 
(multidimensional) Schrodinger pair of canonical observables. Define 



in Hilbert space 



operators pf^ and q^^ 



d, (called 



momentum and position operators, respectively) by the formulae 

{pf^ij){xi,...,Xd) = {Dkilj){xi,...,Xd), 
{qf^'ilj){xi, ...,Xd) = -ixki^ixi, ...,Xd), 
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where denotes the k-i\i partial derivative. The operators p^.^"* and qf"* 
are unbounded densely defined and selfadjoint; moreover, they satisfy 
the commutation relations 

(6) [Pf,pf] = [.f,gf] = 0, [pi^,f] = -.4,l, 

where for operators A, B on L^(]R^) 

[A, B] = AB^ BA, 

and 1 stands for the identity operator (observe that since pf^ and q^^^ 
are densely defined, relations (E]) are assumed to hold only on a dense 
subspace of L'^{R'^)). 

The pair {p^^\ q^^^) = {{Pi'\ Qi^^), ■ ■ ■ , {Pd'\ Qd^^)) is called the Schro- 
dinger pair of canonical observables. Putting 



V''^\xi,yi, Xd, yd) = exp <j z ^ {xkPk^ + VkQ^) 



{xi,yi, ...,Xd,yd) e 



k=l 
n2d 



we easily see that z V^^\z) is a projective unitary representation of 
the group R^'^ on Ti, satisfying the Weyl-Segal commutation relations 
([1]). Now if T^°) is a probability operator on L^(M'^) (we use a super- 
script when referring to the space L^(M'^)) then its characteristic 
function at the point A^z for An given by the formula ([5]) equals to 



f W (A^z) = tr r(°)y W (Ar^z) = tr T^") exp [tj^ 4"^ (^4°^ + ) 

k=l 

= tr T(°) exp {^ f^^ ^(4^^^ ) + l/.(4"^gf )] } 

k=l 

which corresponds to the passing from the multidimensional canonical 
pair {{pf\qf^^), {Pd\qf^)) to the pair 

((a;"'pr,4"'«r).....(«fri°'.«f9r)). 

i.e. each of the component pairs {pfKq^'^) being normed by possibly 
different numbers a^"\ k = 1, . . . ,d. It is worth noting that in the pio- 
neering paper [1] on quantum limit theorems, the central limit theorem 
was formulated just in the language of canonical pairs, though solely 
in the case d = 1 and with the classical scalar norming ai"^ = 
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Coming back to our setup, we have 
lim S4z) = lim e^<^"'^> [f (aS^^xi, aS^^yi, . . . , ai^^x,, y.)]" 

= S{xi,yi, ...,Xd,yd), z = {xi,yi, ...,Xd,yd) G IR^''- 

It was proved in [4J that then the hmit operator S must be Gaussian, i.e. 
iS* is the characteristic function of a Gaussian probabihty distribution on 
M^'^. In the classical commutative situation various sufficient conditions 
on belonging to the domain of attraction of a Gaussian law have been 
obtained — the most celebrated being that of finite variance as in 
the Lindeberg-Levy central limit theorem. It turns out that in the 
quantum case this condition is also necessary. Namely, we shall prove 
the following 

Theorem. Let T be an arbitrary probability operator onTi. T belongs 
to the domain of attraction of a Gaussian probability operator if and 
only if T has finite variance. 

3. Proofs 

We begin with a simple lemma which gives a description of the char- 
acteristic function of Gaussian probability operators in a particular 
case. 

Lemma 1. Let 

f{z) = e-tl'^l'^ 

for some a > 0. f is the characteristic function of some Gaussian 
probability operator if and only if a ^ |. 

Proof. Observe that / is the characteristic function of a Gaussian prob- 
ability measure with the covariance matrix Q = aL From [3", Chapter 
V, $$ 4, 5] (see also [6J) it follows that an arbitrary positive-definite 
2d X 2d matrix Q is the covariance matrix of a Gaussian probability 
operator if and only if the following inequality holds 

(8) {Qz, z) + {Qz', z') ^ A{z, z'), z, z' G M'^ 

which in our case amounts to saying that 

«(lkf + Ik'f) > ^{z,z'), z,z' e R^''. 

The inequality above may be rewritten in the form 

d 

^{axl + ayl + + ay[^ - xuy'k + 2/fc4) ^0, Xfc, x'^, G M. 
k=i 
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It is easily seen that this inequahty holds if and only if for each 
k = 1, . . . ,d and arbitrary Xk, i/k, x'j^, G M we have 

axl + ayl + ax'^^ + ay'^^ - xty'k + Vkx'k ^ 0, 
which, in turn, is equivalent to the positive definiteness of the matrix 



a -i 
a i 
i a 
-| a 

Since the eigenvalues of this matrix are equal to a ± ^ the conclusion 



follows. 



□ 



We also have the following simple property of the covariance matrix 
of a Gaussian probability operator 

Lemma 2. Let Q he the covariance matrix of a Gaussian probability 
operator. Then Q is non-singular 

Proof. Indeed, assume that Qz' = for some z' ^ M^'^. Then for 
each fixed z G M^°' and an arbitrary t G M we have on account of ([8]) 

{Qz, z) = {Qz, z) + {Qitz'), itz')) ^ A(z, tz') = tA{z, z'), 

which is clearly impossible. □ 

The following proposition provides estimation on the coefficients of 
the norming matrices. 

Proposition 3. Let {An} be an admissible sequence of matrices of the 
form ([5]). Then 

a^^"^ ^ for each k = 1, . . . ,d. 



n 



Proof. Let Ti = ■ ■ ■ = T„ = T be Gaussian probability operators with 
the characteristic function 



Then 



n 



k=l 



exp 



fc=i 
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which is a Gaussian probabihty operator with covariance matrix 



ai")^ 



(n)2 



al 











a':'' 



a 



(n)2 



Now the inequahty (IHD takes the form 



n 
2 



Yl "fc"^^ i^k + yl + x'k + y'k) ^ Y^ixkv'k - Vkx'k). 

k=l k=l 

Putting x[ = -xi, y[ = yi, Xk = Uk = x'^. = y'l, = for k = 2, d we 
obtain 

(n)2/ 2 I 2\ \ r) 

which means that the matrix 

.(n)2 



na] 



-1 



-1 na 



(n)2 



is positive definite. Consequently, 



I.e. 



2 (")4 ^ 1 



(n) \ 



By the same token we obtain the required inequahties for k = 2, . . . , d. 

□ 

The next lemma is a known classical result from the theory of do- 
mains of attraction (cf. [2,, Chapter IX, $ 8]). 

Lemma 4. Let u be a probability measure belonging to the domain of 
attraction of a Gaussian law, i.e. there are constants 6^ > 0, c„ G M 
such that 

lim e**"" iHbJ)] " = e^*™- , t eR, 

n~>oo 

for some m G M, cr > 0. Ifbn^^, then v has finite variance, i.e. 



\^ v{d\) < oo. 
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Proof. We shall follow [2\. First, note that the theory of limit laws 
yields 6^ — * 0. Fix an arbitrary x > and denote 

U{x) = I X^u{dX). 

According to formula (8.12) in [2, Chapter IX, $ 8, Theorem la] we 
have 



MB-'- 



for some constant c. (We warn the reader that there is a difference in 
the notation employed in [2] and here, namely, we use 6„ for what in 
[2] is denoted by ^ and c„ for what in [2] is denoted by Since 



nbi ^ 1, 



we get 



X^i^idX) = lim / X^u{dX) = lim u(^) < oo. 

On 

□ 

Now we are in a position to prove our theorem. 

Proof of the Theorem. Necessity. Assume that for some probability op- 
erator T, a sequence {zn} of vectors from M^"^ and a sequence {An} of 
admissible matrices of form ([5]) we have 

lim e^<-'^) [f (aS")xi, aS"Vi, • • • , a^^\,, afy,)]"- 

(9) 

= S{xi,yi, ...,Xd,yd), 
for each z = {xi,yi, ...,Xd,yd) € K^'^. Since S is Gaussian 
(10) S{z) = e*<^0'^)-i<«^'^) 

for some Zq G M^"^ and covariance matrix Q. Let fi^ be the probability 
measure defined by the formula ([2]). Our aim consists in showing that 
fiz has finite second moment. We have 

/oo 
e''^ tiTE,{dX) =fl,{t), teR. 
-oo 

Fix z = {xi,yi, Xd, Vd) G K^'^, and put 

zi = (xi,?/i,0, ... ,0), 22 = (0,0,2;2,?/2,0, ... ,0), ... , 
Zd = (0, . . . ,0,Xd,2/d). 
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Assume for a while that for each k = 1, . . . ,d, Zk 7^ 0. We have on 
account of P). ([TUD and ([TT]) 



hm \%Aa^h)re'"^'"^'>'^ = Mm Ifia^j^kz, 

n— >oo n^oo L 

1+2, 



n _ 



-e 



it{zo,Zi:)- ^t^ {Qzk,Zk) 



From Proposition [3] and Lemma H] we obtain that all the measures 
fizki k = 1, . . . ,d have finite second moments, 

Of course, the same is true if = 0, because then fi^^ is the Dirac 
measure concentrated at zero. 

From the commutation relations ([T]) it follows that the unitary groups 
{V{tz) : t e R}, {V{tiZi) : ti e R},..., {V{tdZd) : U e R} form a 
commuting system of operators; moreover, 

(12) e'*^(") = V{tz) = V{tzi) ■ . . . ■ V{tZd) = e'*^(""^) ■ . . . ■ e''^^^'''^ 

for each t G M. It follows that there is a spectral measure F and Borel 
functions f,fk,k = l,...,d such that 

/oo poo 
/(A)F(dA), R{zk)= / /fc(A)F(dA), 
■00 J —00 

and the equality (fT2|l yields 

/(A) = /i(A) + --- + /,(A). 
Furthermore, substituting t = /(A) we obtain 

/oo poo 
fiX)FidX)= / t(/oF)(c^t), 
■00 J —00 

where 

(/oF)(A) = F(ri(A)), AgS(R). 
On the other hand we have 



/oo 
\E,{dX) 
-00 



and the uniqueness of the spectral decomposition yields the equality 

E, = foF. 
By the same token we obtain the equalities 

E2^= fkoF, k = l,...,d. 
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Consequently, we get 

/oo /»oo 
ii TE,{dt) = I ^rT{f o F){dt) 
■oo J —oo 

/CO 
fix) tiTFidX), 
-oo 

and analogously 

POO 

m2ifi2,)= fix) tr TF{dX), k = l,...,d. 



Finally, we have 

f{X) = [/i(A) + • ■ ■ + /,(A)]^ ^ d[mX) + ... + /J(A)] ^ 
yielding 



/OO poo ^ 

fix) tiTFidX) ^ / dJ2fkW tiTFidX) 
■OO J -OO 

^ poo ^ 

= dJ2 fkW tiTFidX) = dJ^Mf^z,) < oo, 

k=i -^-"^ k=i 

which ends the proof of necessity. 

Sufficiency. A proof of sufficiency is essentially contained in [Ij, 
however, since the setup of [1] is different from the one adopted in our 
work and since some considerations about centring should be taken 
into account we present a short proof. Let T be a probability operator 
having finite variance. Take 

1 



Then the sequence of norming matrices {An} reduce to the sequence 
of numbers {"^}, and from [6,, Proposition 2.5] it follows that this se- 
quence is admissible (this can also be checked straightforwardly, namely, 
it is to be verified that the function 



k=i V 



is A-positive definite for arbitrary probability operators Ti,...,T„). 
For an arbitrary z G M."^"^, let as before fi^ be the probability measure 
defined by the formula ([2]). The mean value of fXz equals to mf (z); 
moreover, it is pointed out in |3], Chapter V, $ 4] that mf is a lin- 
ear function of z, which can be checked using the known formula for 
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moments of a probability measure: 

mf(.) = -4Mt)i„ = = At "^^(*^' 

Consequently, there are vectors Zn £ IR^'' such that 

{zn^ z) — —m\{z)\/n for each z e R^'^. 

We have 



t 



/n J ' \\Jni 

Prom the classical Lindeberg-Levy central limit theorem it follows that 



t 



lim e~**'"i^(^)^ 

n— »oo 



t 



for some cr^ > 0, which means that 



hm e**<^-^> 

n— »oo 

Putting t = 1 we get 

hm e*^^"'^^ 



= e 2^ 



e 2^^^^ . 



and the existence of the limit on the left hand side means that on 
the right hand side we have the characteristic function of a Gaussian 
probability operator which finishes the proof. □ 

References 

[1] CD. Cushen and R.L. Hudson, A quantum mechanic central limit theorem, J. 

Appl. Probab. 8(1971), 454-469. 
[2] W. Feller, An Introduction to Probability Theory vol. II (in Polish), PWN, 

Warszawa, 1978. 

[3] A.S. Holevo Probabilistic and Statistical Aspects of Quantum theory, North- 
Holland, Amsterdam-New York-Oxford, 1982. 

[4] K. Lubnaucr, Operator-stable quantum probability operators, J. Math. Sciences 
106(2001), 2713-2718. 

[5] K. Lubnauer, Quantum limit theorems, Studia Math. 164(2004), 103-120. 

[6] K. Urbanik, Non-commutative probability limit theorems, Studia Math. 
78(1984), 59-75. 

Faculty of Mathematics and Computer Science, Lodz University, ul. 
S. Banacha 22, 90-238 Lodz, Poland 

E-mail address: lubnauer@math.uni.lodz.pl 



E-mail address: anluczak@math.uni.lodz.pl 



